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I. INTRODUCTION 



The g-state Potts antiferromagnet (AF) |],|2| exhibits nonzero ground state entropy, S > 
(without frustration) for sufficiently large q on a given graph or lattice. This is equivalent 
to a ground state degeneracy per site W > 1, since S = k B hiW. Such nonzero ground 
state entropy is important as an exception to the third law of thermodynamics |||4[]. The 
zero-temperature partition function of the above-mentioned g-state Potts antiferromagnet 
on a graph G satisfies 

Z(G,q,T = 0) PAF = P(G,q) (1.1) 

where P(G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color 
HIJ. Thus, formally, 

W({G},q)= limPfG.g) 1 /" (1.2) 
where n = v(G) is the number of vertices of G, and we use the notation 

{G} = lim G (1.3) 



n^oo 



The minimum number of colors needed for this coloring of G is called its chromatic number, 
x{G). At certain special points q s (typically q s — 0, 1, .., x(G)), one has the noncommutativ- 
ity of limits 

lim lim P{G,q) 1/n + lim lim P{G,q) 1/n (1.4) 

q— >q s n— +oo n— >oo q^q s 

and hence it is necessary to specify the order of the limits in the definition of ^({G}, q s ) J7|. 
As in [0], we shall use the first order of limits here; this has the advantage of removing certain 
isolated discontinuities that could be present in the W with the opposite order of limits. In 
addition to Refs. @- 0, some other previous related works include Refs. ||- [piR . Since 
P(G, q) is a polynomial, one can generalize q from Z + to 1R and indeed to C. ^({G}, q) is a 
real analytic function for real q down to a minimum value, q c ({G}) fl7lJl4|. For a given {G}, 



we denote the continuous locus of non-analyticities of W as B. This locus B forms as the 
accumulation set of the zeros of P(G, q) (called the chromatic zeros of G) as n — ► oo and 
satisfies B(q) = B(q*). A fundamental question concerning the Potts antiferromagnet is the 
form of this locus for a given graph family or lattice and, in particular, the maximal region 
in the complex q plane to which one can analytically continue the function iy({G}, q) from 
physical values where there is nonzero ground state entropy, i.e., W > 1. We denote this 
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region as R\. Further, we denote as q c the maximal point where B intersects the real axis, 
which can occur via B crossing this axis or via a line segment of B lying along this axis. 
Although W({G},q) has been calculated exactly for the triangular lattice [|ll|] and a 



number of different families of graphs @, pTSfl - p3| , it has never been calculated exactly (as 



a general function of q) for any other 2D (or higher-dimensional) lattice. Accordingly, one 
useful procedure is to calculate W exactly on infinite-length, finite-width strips of various 
2D lattices. As the width increases, the values of the resulting W functions approach those 
of the W functions on the respective 2D lattices. We have carried out this procedure for 
strips with free boundary conditions in the longitudinal direction (in which the length goes to 
infinity) and free and periodic boundary conditions in the transverse direction |T9| , |22[| . For a 
given lattice, the calculation of P and W ona. cyclic strip of a lattice, i.e a strip with periodic 
boundary conditions in the longitudinal direction (and some choice of boundary conditions 
in the transverse direction) is more difficult and the results more complicated than is the 
case for the corresponding strip with free longitudinal boundary conditions (and the same 
transverse boundary conditions). Despite the greater difficulty, the calculations of W on 
cyclic strip graphs are worth the effort, because, as explained further below, the resultant 
W functions exhibit, for finite L y , certain analytic features expected for the W functions on 
the infinite 2D lattice, whereas for strips with free longitudinal boundary conditions, these 
features are absent (approached in the limit L y — > oo). 

Accordingly, we have carried out exact calculations of the chromatic polynomials P and, 
in the infinite-length limit, the resultant functions W, for cyclic strips of the square and 
kagome lattices, and we report the results here.Q Since the W functions that we have 
computed for these strips exhibit the analytic properties expected for the W functions on 
the respective 2D lattices, they constitute, in this sense, the closest exact results that one 
has to these W functions. We also include results for some other families of cyclic strips. 
For all of these strips the boundary conditions in the transverse direction are free. 

From our earlier exact calculations of If on a number of families of graphs we have 
inferred several general results on B: (i) for a graph G with well-defined lattice structure, a 
sufficient condition for B to separate the q plane into different regions is that G contains at 
least one global circuit, defined as a route following a lattice direction which has the topology 
of S 1 and a length l g c _ that goes to infinity as n — ► oo |19| .[] 



1 A brief report of some of the results for W is given in Ref. p5[ . 

2 Some families of graphs that do not have regular lattice directions have noncompact loci B that separate 
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For a d- dimensional lattice graph, the existence of global circuits is equivalent to having 
periodic boundary conditions (BC's) in at least one direction. Further, (ii) the general 
condition for a family {G} to have a locus B that is noncompact (unbounded) in the q plane 
16fl shows that a sufficient (not necessary) condition for {G} to have a compact, bounded 



locus B is that it is a regular lattice . The third and fourth general features are that 

for graphs that (a) contain global circuits, (b) cannot be written in the form G = K p + H 
jlSfPi an d (c) have compact B, we have observed that B (iii) passes through q = and (iv) 
crosses the positive real axis, thereby always defining a q c . As noted, the advantage of cyclic 
strip graphs is that these properties are present for each finite L y rather than only being 
approached in the limit L y — > oo as for open strips. 



Our present results also strengthen the evidence for our earlier conjectures ||19| , |23|| that 
on a graph with well-defined lattice structure, a necessary property for there to be chromatic 
zeros and, in the n — > oo limit, a locus B including support for Re(q) < 0, is that the graph 
has at least one global circuit. (This is known not to be a sufficient property, as shown, e.g., 
by the circuit and ladder graphs, which have such global circuits, but whose chromatic zeros 
and loci B have support only for Re(q) > 0.) 

II. GENERAL PROPERTIES AND CALCULATIONAL METHOD 

A general form for the chromatic polynomial of an n-vertex strip graph G of the type 
considered here (and more generally, a graph built up recursively by the successive addition 
of a given subgraph) is 

N a 

P{G,q) = Y, c M){aM)) tin (2-1) 



the q plane into different regions |p^| , |l8| , p3| . 

3 The complete graph on p vertices, denoted K pi is the graph in which every vertex is adjacent to every 
other vertex. The "join" of graphs G\ and G*2, denoted G\ + G2, is defined by adding bonds linking each 
vertex of G\ to each vertex in G2. Graph families with B not including q = are given in Jl5|,^6|fl^,|2^] . 
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where Cj(q) and aj(q) are certain functions of q. Here the aj(q) and Cj(q) are independent 
of n, and tj is a G-dependent constant. As before 0, we define a term a^(g) as "leading" 
(£) if it dominates the n — > oo limit of P(G,q); in particular, if iV a > 2, then it satisfies 
|a^(g)| > for j ^ £, so that \W\ = \ae\ te . The locus B occurs where there is a 

nonanalytic change in W as the leading terms ai in eq. (|2.1|) changes. Note that a term 
of the form a\ — ±1 may be absent. For some families of graphs, the Cj(q) and cij(q) are 
polynomials. However, there are also many families of graphs for which Cj(q) and a,j(q) are 
not polynomial, but instead, are algebraic, functions of q; for these families the property that 
the chromatic polynomial is, in fact, a polynomial of q, is not manifest from the expression 

For a strip graph G s our calculational method is to use the deletion-contraction theorem 
H iteratively to obtain a set of linear equations which we then solve, to get a generating 
function T(G s ,q, x); the chromatic polynomials for the strip of length L x = m are then 
determined as the coefficients in a Taylor series expansion of this generating function in an 
auxiliary variable x about x = 0: 

oo 

T(G s ,q,x)= P((G s )m,q)x m - mo (2.2) 

m=mo 

where m depends on the type of strip graph G s and is naturally chosen as the minimal value 
of m for which the graph is well defined. The generating functions T(G S , q, x) are rational 
functions of the form 

r(G„,,x) = ^4 (2.3) 
V(G s ,q,x) 



with 



N{G s ,q,x) = Y< A GsMx j (2.4) 

j=0 



and 



The labelling convention for the summation index j 
differs from the convention we used in Ref. [Q, where 
Ci(a\) tin was labelled simply as cq since cq was either 
with n like {const.) 71 with |consi.| > 1. 



in eq. ( |2.1| ) is the same as that in our Ref. this 
j started at 0, and where the term denoted here as 
1 or proportional to (—1)™, and hence did not grow 
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T>(G s ,q,x) = l + ^2b GsJ (q)xi (2.5) 

3=1 

where the Aq s ^ and bc a ,i are polynomials in q (with no common factors) fj and the degrees 
of the numerator and denominator, as polynomials in the auxiliary variable x, are 

djv = deg x {N) (2.6) 

and 

d v = deg x {V) (2.7) 
Writing the denominator of the generating function T(G s ,q,x) in factorized form, we have 

V(G s ,q,x) = H(l-\ GsJ (q)x) (2.8) 
j=i 

(c.f. eq. (4.2) in Ref. [|19[ or (2.9) in Ref. [^H). The chromatic polynomials can equivalently 
be calculated from a corresponding recurrence relation (see eq. (2.25) in Ref. pTf). 
We define the polynomial 

Dtiq) = = ^d-^ = b-^'C ~. i y~ 2 " <2 ' 9) 

where a = q — 1 and P(Ck, q) is the well-known chromatic polynomial for the circuit (cyclic) 
graph Ck with k vertices, 

P(C k ,q) = a k + (-l) k a (2.10) 

In Table 1 we list the values of dj^ and dv for various families of lattice strip graphs 
considered in this paper and in our previous ones. The form of Af and value of d^f, (but 
not T> or d-p) depend on one's convention for m in eq. ( p. 2D ; the values given in Table 1 
correspond to the conventions that we use in this paper and our previous ones. 



In eq. (2.17) of Ref. |2T] we gave a general formula that enables one to calculate the 



chromatic polynomial in the closed form (|2.1J ) from the generating function or vice versa. 



5 As in Ref. |2l)|, we use the notation Aq s i rather than our earlier notation ag Jit in Ref. pil in order to 



avoid confusion with the terms a,j in the general formula (2.1). 
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Since the denominators of generating functions for wider strips of various lattices G s yield 
Ag S) /s in eq. (|2.8|) that are algebraic, but not polynomial, functions of q, it can be more 
convenient to obtain the chromatic polynomials directly from the expansion of the generating 
functions. However, for the simplest types of strips, the Ag sJ 's are polynomial, and in these 
cases the closed-form expression ( |2.1|) is the most convenient. In the case of cyclic strips 
of the square lattice, it is worthwhile to give the generating functions for the simpler cases 
L y — 1 (the circuit graph) and L y = 2 (the cyclic ladder graph) as a comparison with the 
much more complicated results for the L y = 3 case that we are reporting here. For these 
cyclic strips one can write 

oo 

r(sq(L y ,cyc.),q,x) y ,cyc.) m ,q)x m 2 (2-11) 

m=2 

As indicated, it is natural to take mo = 2 for the cyclic strip of the square lattice (and for 
the cyclic strip of the triangular lattice to be discussed later) because the lowest value of m 
is m = 2, for which the respective strips reduce to a single column of squares or triangles. 
For the minimal case, L y = 1, the generating function that yields, as the x m ~ 2 coefficient, 
the chromatic polynomial P(C m ,q), is 

T(sq(L y = 1, eye), q, x) = f^'^L, , ( 2 - 12 ) 

[l—(q — l)x\(l + x) 

As an example of the mo-dependence of M , we note that if one used a definition with m = 0, 
i.e., T(sq(L y ,cyc.),q,x) = J2^ =0 P(sq(L y ,cyc.) m ,q)x m , (where the m = 0, 1 coefficients are 
only formal) then 

fV (T 1 ^ ^ q[l-(q-2)x] 1 , g-1 

T(sq(L y = 1, eye), q, x) = {1 _ iq _ 1)x]{1+x) = 1 _ {q _ 1)x + — x 

oo 

= l) m + {q- l)(-l) m x m (2.13) 

m=0 

which, of course, agrees with the expression for P(sq(L y = l) m , q) = P(C m , q) given in eq. 
( p.lO| ). For L y = 2, from the known chromatic polynomial |J, we determine the correspond- 
ing generating function for eq. ([2.1 1|) (with the convention mo = 2) to be 

T(sa(L -2 eve) ax)- ^ ~ - 2(g - l)x - (g - l)(g - 3)D 4 x 2 ] 

1 {sq(L y - 2, eye), q, x) - ^ _ _ x)[1 _ {1 _ q)x][1 _ (3 _ q)x] ^) 

where = q 2 — 3q+3. As will be seen from our new results, the complexity of the generating 
function increases considerably when one goes from width L y = 2 to L y = 3. 

TABLE I. Properties of generating functions for strip graphs G s . 
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For a given G S} in region Ri, the function W is given by 

W = (XmaxY (2.15) 

For q in regions other than Ri, one can only determine the magnitude q)\ unam- 

biguously 0; this is given by 

\W\ = |A ma ,|* (2.16) 

where A max denotes the A with maximal magnitude in the respective region and 

m 

t=— 2.17 

n 

For the square and triangular strips, t = 1/L y while for the kagome strip considered here, 
t = 1/5. Thus, W and B are determined by V, independent of Af. Of course, the de- 
termination of TV is necessary to calculate the actual chromatic polynomials for particular 
finite-length strips. Besides their intrinsic interest in the context of mathematical graph the- 
ory, these chromatic polynomials for finite strips are important because from them one can 
calculate their zeros and thereby check how rapidly these approach the asymptotic n — > oo 
accumulation set B. (Indeed, for some other families of graphs, we have shown that the 



chromatic zeros lie precisely on the asymptotic loci even for finite n ||15|| .) 

The organization of the paper is as follows. Our calculations for the strips of the square, 
kagome, and triangular lattices are given in sections 3-5 (see also the appendix). In section 
6 we include results for certain homeomorphic expansions of square strips. Our concluding 
remarks are given in section 7. 



III. CYCLIC STRIP OF THE SQUARE LATTICE WITH L y = 3 

We consider the cyclic strip of the square (sq) lattice of length L x = m vertices and width 
L y = 3 vertices. This has n = L x L y = 3m vertices, e = 5m edges (for the nondegenerate 
cases m > 2), and chromatic number (for any L y ) 

{2 for even m 
(3.1) 
3 for odd m 

(For the minimal case m = 2, the periodic longitudinal boundary conditions lead the strip 
to degenerate to a single tower of two squares.) We use a generating function of the form 

oo 

T(sq(L y = 3, eye), q, x) = ^ P ( S( l( L y = 3 > cyc.) m , q)x m ~ 2 (3.2) 

m=2 



s 



where P(sq(L y = 3, cyc.) m , q) is the chromatic polynomial for the strip of length m vertices. 
We obtain a generating function of the form (2J3) with (]2.4j) and (|2.5|) , where deg x (M) = 
8 and deg x (T>) = 10. Clearly this is considerably more complicated than the generating 
function for the strip with the same width, L y = 3, but with free longitudinal boundary 
conditions (and again free transverse boundary conditions), which had deg x (M) = 1 and 
deg x (T>) = 2. This and other comparisons are shown in Table 1. 
We have calculated the denominator to be 



V(sq(L y = 3), q, x) = [1 + x] [1 - (q - l)x] [1 - (q - 2)x] [1 - {q - A)x\ [1 + (q - 2) 2 x] 



where 



x [1 + b aqjll x + b sqtl2 x 2 ][l + b sqt 2ix + b sqt22 x 2 + b sg:23 x 3 



(3.3) 



J sq,U 



-(g-2)(g 2 -3g + 5) 



(3.4) 



b sqA2 = (q-l)(q 3 -6q 2 + 13q-ll) 



b sq ,2i = 2q 2 - 9q + 12 



(3.5) 
(3.6) 



b sq ,22 = q 4 ~ 10g 3 + 36q 2 - 56q + 31 



(3.7) 



Thus, 



b sq>2 3 = -{q - 1)(? 4 - 9g 3 + 29g 2 - 40g + 22) 
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V(sq(L y = 3, eye), q, x) = JJ(1 - X sqd 

3=1 



x 



(3.1 



(3.9) 



where the first five A's can be read off immediately from the linear factors in eq. (|3.3| ) 

1 



and X S qj, j = 8, 9, 10 are the roots of the cubic 

£ 3 + bgq^i^ 2 + 6 S g 5 22^ + frsg,23 







(3.10) 



(3.11) 



From this we have calculated the locus B, shown in Fig. |l]. This locus separates the q 
plane into seven regions: R\ including the real interval q > q c ; R 2 , to the left of Ri, R3 
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including the real interval < q < 2; two complex-conjugate (c.c.) phases R±,R\ centered 
at q ~ 2.4 ± 0.92; and two additional quite narrow, sliver-like c.c. phases R§, R^ centered at 
q ~ 1.95 ± 1.75z. In region R±, the dominant A is A S(?i 6 = \ a q,i,+, so that 



W(sq(L y = 3, eye), q) = 2" 1 / 3 [(g - 2)(g 2 - 3g + 5) + 

(g 2 - 5g + 7)(g 4 - 5g 3 + llg 2 - 12g + 
In region R 2 , the dominant A is q — 4, so 

|W(sg(L y = 3, eye.), g)| = |?- 4 l 1/3 for g e R 2 



1/2 



1/3 



for q E R\ (3.12) 



(3.13) 



The value of q c for this strip is given as the real solution to the equation {Wr^ = \Wr 2 \, 
which reduces to 



yielding 



2g 4 - 16g 3 + 51g 2 - 86g + 67 = 



q c (sq(L y = 3, eye.)) = 2.33654 



(3.14) 



(3.15) 



At this point W = 1.18487. In region R%, the dominant A is the maximal root of the 
cubic ( p. 11| ). The degeneracy condition \Wr 2 \ = \Wr 3 \ yields the value q = 2 for the point 
separating R 3 and R 2 on the real axis. Note that the portion of the boundary B separating 
regions R 2 and R% is an arc of the circle |g — 4| = 2. In the c.c. pairs of regions (R4, R£) and 



(i?5,i?g) the dominant A's are the other two roots of the cubic ( 3.11 ). As discussed before, 
these results depend only on D. As is evident in Fig. [1], B contains several multiple points, 
forming conjugate complex pairs. The property that there is no multiple point at q c follows 
because only two of the A's become degenerate in magnitude at this point. This may be 
contrasted with the situation for the L y = 2 cyclic strip, where three leading terms become 
degenerate, leading to a multiple point, at q c (sq(L y = 2, eye.)) = 2, at which four curves 
come together |7[|. As noted above, these results depend only on V. 

The coefficient functions for the numerator are rather lengthy and hence are given in the 
Appendix. Together with the above results for T>, these yield, via eq. ( |3.2| ), the chromatic 
polynomials for the cyclic strip of the square lattice of arbitrary length. From the generating 
function, using eq. (2.17) of Ref. |2l[] (with m — > m — 1 to match our current notational 
conventions), we calculate the chromatic polynomial 



P{sq{L y = 3, cyc.) m , q) = (g 3 - 5g 2 + 6g - 1)(-1) 
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+ (g 2 -3g + l) (g-l)- + ( g -2) m + (g-4) m + (q - l)[-(q - 2) 2 } m + (X sqfi ) m + (X sqJ ) 



+(?-!) (\s g ,s) m + (K q , 9 ) m + (A S9)10 ) 7 



This can also be derived directly from T> without using M by matching with calculations of 

or eq. (|27l 

5g 2 + 6g — 1 



P for the first few values of m. Thus, in the notation of eq. (2.17) of Ref. [^T 

3 



here, a sqJ 



A 



sq,ji Csq,l 



q 



-> c sq,j 



q — 3q + 1 for j 



2; 3, 4, C S q j 



(1 



for j = 5, 8, 9, 10, and c sq j = 1 for j = 6, 7. In Fig. [I] we show chromatic zeros for a strip of 
length m = 20 and thus n = 60 vertices, for comparison with the n — > oo accumulation set 
£>. Besides the fact that these zeros lie close to £>, we observe that the density is largest on 
the right-hand side of the "main" curve in B and is rather small for the curve bounding the 
region R2 on the left. 



(3.16) 



2.0 



1.0 



0.0 



-1.0 



-2.0 1 1 1 ■ 1 ■ 1 1 1 1 1 1 

-1.0 0.0 1.0 2.0 3.0 4.0 

Re(q) 

FIG. 1. Boundary B in the q plane for W function for a cyclic strip of the square lattice with width 
L y = 3. Chromatic zeros for m = 20 and thus n = 60 vertices are shown for comparison. 

An interesting feature of this family of graphs is that the chromatic zeros and their 
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accumulation set B have support for Re(q) < 0. Chromatic zeros with negative real parts 
occur first for m = 10, i.e., n = 30. 



IV. CYCLIC STRIP OF THE KAGOME LATTICE 

We next consider a cyclic strip of the kagome lattice comprised of m hexagons with each 
pair sharing two triangles as adjacent polygons (as shown for the open strip in Fig. 1 in Ref. 
19fl). This has n = 5m vertices and e = 8m edges. The chromatic number is \ = 3. We use 



a generating function of the form 

oo 

T(kg(L y = 2,cyc.),q,x) = ^ P(kg(cyc.) m , q)x m ~ 1 (4.1) 



m=l 



where P(kg(cyc.) m , q) is the chromatic polynomial for the cyclic strip with m hexagons. The 
generating function V has deg x (N) = 8 and deg x (T>) = 9, and is thus much more complicated 
than the generating function for the kagome strip of the same width with free longitudinal 
boundary conditions, which had deg x (Af) = 1 and deg x {V) = 2 (see Table 1). We have 
calculated 

V(kg(L y ,cyc.),q,x) = (1 + h g>n x + b kgtl2 x 2 )(l + b kg ^ix + b kg ^2X 2 ){l + b kgi31 x + b kgfi2 x 2 )x 

[1 - (q - 2)x][l -{q- A)x\[l - (q - l)(q - 2) 2 x] (4.2) 



where 



bka.li = -(<? " 2)(g 4 - 6g 3 + 14g 2 - 16q + 10) (4.3) 

b kg , 12 = (q-lf(q-2) 3 (4.4) 

b kg ,2i = ~q 3 + 7q 2 - Wq + 20 (4.5) 

b k9 ,22 = (g-l)(g-2) 3 (4.6) 

b kg:31 = ll-9q + 2q 2 (4.7) 



b kg<32 = -(q-l)(q-2) 2 
12 



(4.8) 



We thus find, for 1 < j < 6, 



^kg,j — ^kg,k,± 



1 



-h g ,ki ± (&L,fci - 46 fe9ifc2 ) 1/2 



Afc g ,7 — Q — 2 



V,s = q - 4 



A fe9 , 9 = (g-l)(g-2) 5 



, fc=l,2,3 



(4.9) 
(4.10) 
(4.11) 
(4.12) 



2.0 



1.0 



-1.0 




-2.0 1 1 1 1 1 1 1 1 1 1 1 1 

-1.0 0.0 1.0 2.0 3.0 4.0 

Re(q) 

FIG. 2. Locus B for W for oo x 2 cyclic strip of the kagome lattice. Chromatic zeros for m = 10 (n = 50) 
are also shown for comparison. 
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The coefficient functions for the numerator are rather lengthy and hence are given in 
the Appendix. From the generating function, via eq. (2.17) of Ref. |21|, or directly from T> 
by matching results with computations of P for the first few values of m, we calculate the 
chromatic polynomial 



P(kg(L y = 2) m , q) = {q 2 -3q + l)(q- A) m + (q - 1) (q - l)(q - 2 



+ 



(A fc9 ,i) m + (\ kg , 2 ) m +(<?-!) (A fc9 , 3 ) m + (A fcffi4 ) 



(4.13) 



There are three A's in T> that do not appear in P, namely, q — 2 and Xkg,(5,6) = ^kg,3,±- in 
the context of the general formula, eq. (2.17) of Ref. ||21|| , this can be seen as a result of 
the vanishing of the corresponding coefficient functions Cj in eq. (2.19) of that work. Fig. ||] 
shows the boundary B in the q plane for the W function. Again, B divides the q plane into 
several regions (Fig. 2). In region Ri, W is determined by Afc S) i = Xk g ,i,+- 



W({G kg{Ly=2) },q) = 2- 1 /\q-2) 1 



/5 



q A - 6q 3 + Uq 2 - lQq + 10 



+ 



q 8 - 12g 7 + 64<f - 200<f + 404g 4 - 548g 3 + 500g 2 - 292g + 92 



1/2 



1/5 



(4.14) 



As q decreases through q c , \W\ switches to \q— 4I 1 / 5 ; solving the degeneracy equation |Afc Si i| = 
|4 - q\ yields q c = 2.62421 (at which point W = 1.06588). This is within about 10 % of the 
inferred exact value q c = 3 for the 2D kagome lattice ||22|| . It is impressive that an infinite 
strip of width L y = 2 yields a q c this close to the value for the full 2D lattice. In the regions 
including the intervals 2 < q < q c and < q < 2, \W\ — \q — 4| 1//5 and \W\ = |A fc9) 4| 1//5 , 
respectively. It is interesting to observe that the boundary B includes a part with Re(q) < 0. 
Of the nine coefficient functions in the polynomial in x in V, the two that are equal to the 
coefficient functions for the open strip of the kagome lattice of the same width are &%,n 
and bkg,n- Thus, the root Xk g ,i is leading in region R\ for both open and cyclic longitudinal 
boundary conditions. 

To compare with the asymptotic locus B, we have calculated the zeros of the chromatic 
polynomial for a finite strip with m = 10, i.e., n = 50 vertices; these are displayed in Fig. 
(|). As in the case of the L y = 3 cyclic strip of the square graph, this cyclic strip of the 
kagome lattice has the property that the chromatic zeros and their accumulation set B have 
support for Re(q) < 0. Chromatic zeros with negative real parts occur first for m = 5, i.e., 
n = 25). 
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V. CYCLIC STRIP OF THE TRIANGULAR LATTICE WITH L y = 2 



In this section we consider the cyclic strip of the triangular (t) lattice of length L x = 
m vertices and width L y = 2 vertices. This has n = L x L y = 2m vertices and, for the 
nondegenerate cases m > 3, e = 4m edges (bonds). We consider both the cases of periodic 
and twisted periodic longitudinal boundary conditions, and we denote the correponding 
strips as cyclic and twisted cyclic (= Mobius, i.e., opposite ends are identified with reversed 
orientation). For m = 2 the cyclic and twisted cyclic strips both degenerate to the complete 
graph on 4 vertices, K±, with e = 6 edges. The cyclic strip has chromatic number 

{3 if m — mod 3 
(5.1) 
4 if m — 1 orm = 2 mod 3 

whereas the twisted cyclic strip has 

x{t(L y = 2, L x — m), twisted eye.) = 4 (5.2) 
For these cyclic strips we use generating functions of the form 

oo 

T(t(L y = 2, eye), q,x) = J2 P (t(L y , cyc.) m , q)x m - 2 (5.3) 

m=2 

oo 

r(t(L y = 2, twisted cyc.),q,x) = P(t(L y , twisted cyc.) m ,q)x m ~ 2 (5.4) 

m=2 

where P(t(L y , cyc.) m , q) and P(t(L y , twisted cyc.) m , q) are the chromatic polynomials for the 
respective strips. For these strips, we find that deg x {T>) = 4 (see Table 1). The degrees in x 
of the numerator for the cyclic and twisted cyclic strips are, respectively, deg x {H) = 2 and 
deg x (J\f) = 0. We calculate 

b t ,i = ~q 2 + 6g - 10 (5.5) 
^,2 = -(g-3)(2g 2 -9g + ll) (5.6) 
^ i3 = -(g-2)V-6g + 10) (5.7) 
h,4 = (q- 2) 4 (5.8) 
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The denominator of the generating function has the factorized form 

V(t(L y = 2, eye), q, x) = (1 - x) [l - (q - 2) 2 xj [l + (2q - 5)x + (q - 2) V] 

4 

= ri( i - x ^ x ) 

3=1 



where 



A t) i = 1 
A il2 = (g-2) 2 



A t)3 - 2 

A*,4 = 7T 



5 - 2g + v/9-4g ] 
5 _ 2g - v/9-4g ] 



(5.9) 

(5.10) 
(5.11) 

(5.12) 
(5.13) 



2.0 



1.0 



f 0.0 



-1.0 



-2.0 1 — 
-1.0 




0.0 1.0 2.0 

Re(q) 



3.0 



4.0 



FIG. 3. Locus B for the function W for a cyclic strip of the triangular lattice with width L y = 2. 
Chromatic zeros for m = 20 (n = 40 vertices) are shown for comparison. 
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Because the locus B depends only on V and is independent of Af, it is the same for the 
cyclic and twisted cyclic strips; hence, in this discussion, we shall group both together and 
drop the "twisted" label. The locus B is shown in Fig. ||] and separates the complex q plane 
into three different regions. We find that 

qML y = 2,cyc.)) = 3 (5.14) 

This may be compared with the value of q c for the full 2D lattice, which is q c (tri) = 4 [nj ,7)- 
The region R\ contains the real interval q > 3. The other two regions, R 2 and R3, include 
the respective real intervals 2 < q < 3 and < q < 2. We find that the leading A's in regions 
Ri, i?2 and R3 are respectively X ma x = A2, Ai, and A3. It follows that 

W{t(L y = 2, eye), q) = D 3 = q - 2 for q 6 R x (5.15) 



\W(t(L y = 2,eyc.),q)\ = l for q e R 2 



(5.16) 



\W(t(L y = 2,cyc.),q)\ 



5-2q+ v/9-4g 



1/2 



for q G R 3 



(5.17) 



As an algebraic curve, the boundary B has a multiple point of index 2 at q = q c in the 
technical terminology of algebraic geometry [2?| (i.e. two branches of the curve cross at this 
point with two separate tangents). As is clear in Fig. (|3|), for this strip, B has no support for 
Re(q) < 0. This is in accord with our earlier inferences noted above, since the global-circuit 
condition is necessary, but not sufficient, for B to have support for Re(q) < 0. 
For the numerator of the generating function, we find, for the cyclic strip, 



A tfl = P{K Al q) = q{q -!)(<?- 2)(q - 3) 



(5.1* 



A ttl = q(q-l)(q-2f 



(5.19) 



A t , 2 = -q(q - l)(q - 2Y 



(5.20) 



where P(K P , q) is the chromatic polynomial for the complete graph on p vertices, 

P-i 

P(K p , q) = H(q - s) 

s=0 



(5.21) 



For the twisted cyclic strip (for which deg x (Af) = 0) we find 
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A u ,o = A t>0 = P(K 4 , q) 
By the same methods as described above, we obtain 

P(t{L y = 2, cyc.) m , q) = q 2 -3q+l + (A t , 2 ) m + (q - 1) [(A t , 3 ) m + (A M ) 

[m/2] , . 

= q 2 - 3q + 1 + (q - 2) 2m + 2 1 - m {q -1) (5- 2g) m - 2s (9 - 4g) s 

s=0 ^ S ^ 



(5.22) 



(5.23) 



P(t(L y = 2, twisted cyc.) m ,q) = -1 + (At )2 )' 

[(m-l)/2] 



(?-!)(? -3) 



(\ 3 ) m - (A M ) ? 



L V // J / \ 

-l + (g-2) 2m -2 1 -(g-l)(g-3) g ^ J 



(5 - 2g) 



m— 2s— 1 



(9-4g) s (5.24) 



where in the upper limits on the sums, [x] denotes the integral part of x. 



VI. ASYMMETRIC HOMEOMORPHIC EXPANSIONS OF THE SQUARE STRIP 

It is also of interest to study homeomorphic expansions of the cyclic square strip In 
Refs. we analyzed such homeomorphic expansions of a square strip with free longitudinal 
boundary conditions (and free transverse boundary conditions). In particular, let us consider 
the case in which for a L y = 2 square strip, we add ki — 2 vertices to each of the upper 
horizontal bonds and k^ — l vertices to each of the lower horizontal bonds. For k\ = /c 2 , the 
resulting A's are polynomial, and the chromatic polynomials can be written in a relatively 
simple form P^j . Here we consider the case k\ ^ k 2 , which is more complicated. We first 
comment on a wider set of graphs, defined such that for each square in the original strip one 
has a choice as to whether to assign the ki vertices to the upper longitudinal side and k 2 
vertices to the lower side, or vice versa. We define a parity vector 



(01, 2 , •••0n 



(6.1) 



in which the &j = + if for the j'th p-gon, the assignment is (upper, lower) = (ki,k 2 ) and 
0j = — if the assignment is (upper, lower) = (k 2 , k\). We shall denote this strip graph as 



A homeomorphic expansion of a graph consists in the addition of degree-2 vertices to bonds of that graph. 
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{Ch)k lt k 2 ,T,,cyc.,m — H EL ^2^2-2 (G sq(L y =2),cyc.,rnj (6-2) 

It will be convenient to define the special E vector 

E +>m = (+,+, ...,+) (6.3) 

where the dimension m of E can be implicit. Obviously, the (Ch) kl>k ^T.,cyc.,m graphs are 
invariant under a reflection about the longitudinal axis, which in the ki 7^ k% case amounts 
to a simultaneous reversal of the signs of all of the cr/s. Let us define this parity operation 
on E as 

P(E) = -E (6.4) 

Then 

(Ch) k (Ch) (Ch) 

k2,ki,T,,cyc.,m 

(6.5) 

The total number of vertices is 



v((Ch) kl ,k2^,cyc.,m) =(p- 2)m (6.6) 

independent of E, where 

p = h + k 2 (6.7) 

For the strip with free longitudinal boundary conditions, we showed in Ref. that the 
chromatic polynomial P((Ch) kl)k2t ^ m , <?) was (i) independent of E and (ii) depended only 
on the homeomorphic expansion indices k\ and k% through their sum, p. However, these 
properties do not hold for the cyclic strip (C7i)fc 1) fc 2) £, C j/c.,m- This can be understood as a 
consequence of the fact that tying together the ends via the periodic longitudinal boundary 
conditions increases the constraints on the coloring of the graphs. 

To study asymmetric homeomorphic expansions, we shall consider the case (ki,k 2 ) = 
(k, 2) (whence p = k + 2) with E = S +>m , i.e., the family 

(C/i)fc 1 =fc,fc 2 =2,S + , m ,cyc.,m (6.8) 

For this calculation we use a generating function of the form 

00 

T({Ch) k ^ + ,cyc, q,x) = J2 X m ~ 2 P((Ch) k ,2,X + , m ,cyc.,m, ?) (6.9) 

m=2 
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From eq. (|6.6| ), the total number of vertices is n = km. We obtain a generating function of 
the form given by eq. (2J3) with ( |2.4|) and ( |2.5|) , where deg x (N) = 3 and deg x (T>) = 5. The 
denominator of the generating function is 



V = [1 - (-l) fc x](l - D k+2 x)[l + (-l) fc (g - 2)x]x 



1 + 



We thus have 



D h+ i + {-!)%- 2) x+ D k+1 (-l) k (q-2)-(-l) k D k 



(6.10) 



(6.11) 



A 



h,2 



D 



k+2 



(6.12) 



-(-l) fc (g-2) 



(6.13) 



Ah, (4,5) 



i - (-l) fc (g - 2) ± ([D fc+1 - (-l) fc (g - 2)] 2 + 4,D fc (-l) fc ) 1 /2 ( 6 .14) 

(where the subscript h means homeomorphic.) For the chromatic polynomial we obtain 

P((Ch) kt2 ,v +intCyc ., n , q) = (q 2 -3q+ l)(-l) km + (A,, 2 ) m + (q - 1) [(A,, 4 ) m + (A,, 5 ) m ] (6.15) 

Note that A^^ does not occur in P. For illustration, the coefficient functions entering in the 
numerator M are given in the Appendix for k = 3, 4. 

Fig. |] shows the analytic structure of the W function for a cyclic strip of the square 
lattice with homeomorphic expansion parameters (k±, k 2 ) = (a) (2,3) (b) (2,4). In these cases 
we find the respective values (a) q c ~ 2.455 and (b) q c = 2. In the region Ri the leading 
term is \h,2- In these cases the boundaries B do not have any part with Re(q) < 0. 
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VII. CONCLUSIONS 



In this paper, we have given exact calculations of the zero-temperature partition function 
of the g-state Potts antiferromagnet strips of the square, kagome and triangular lattices 
with periodic longitudinal boundary conditions (and free transverse boundary conditions). 
These yield, in the limit of infinite length, the ground-state entropy Sq = ks^aW via eq. 
(L2|). These results are of fundamental interest since this model exhibits nonzero ground 
state entropy Sq > for sufficiently large q and hence is an exception to the third law of 
thermodynamics. We also include results for homeomorphic expansions of the square lattice 
strip. The analytic properties of W(q) are determined and related to zeros of the chromatic 
polynomial for long finite strips. 



This research was supported in part by the NSF grant PHY-97-22101. 
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VIII. APPENDIX 

We list in this appendix the coefficient functions for the numerators of the generating 
functions for the cyclic strips of the square and kagome lattices and for asymmetric homeo- 
morphic expansions of the cyclic square strip. 



A. Strip of the Square Lattice 

For the cyclic strip of the square lattice with L y = 3 we calculate 

A sqfi = q(q-l)(q 2 -3q + 3) 2 (8.1) 

A sqA = -2q(q - l)(g 5 - 6g 4 + llg 3 - 21q + 16) (8.2) 

A sq>2 = -q(q - l)(3g 8 - 45g 7 + 306g 6 - 1242g 5 + 3301g 4 - 5854g 3 

+6703g 2 - 4482g + 1318) (8.3) 

A sqt3 = -2q(q - l)(g 10 - 24g 9 + 250g 8 - 1516g 7 + 5992g 6 - 16230g 5 

+30580g 4 - 39552g 3 + 33530g 2 - 16764g + 3732) (8.4) 

A sq ,4 = q(q - l)(7g n - 159g 10 + 1636g 9 - 10077g 8 + 41328g 7 - 118548g 6 

+242569g 5 - 353514g 4 + 358776g 3 - 240824g 2 + 95932g - 17119) (8.5) 

A sq ^ = -2q(q - l) 2 (4g n - 91g 10 + 935g 9 - 5728g 8 + 23249g 7 - 65620g 6 

+ 131260g 5 - 185622g 4 + 181164g 3 - 115604g 2 + 43126g - 7070) (8.6) 

Asqfi = q(q ~ l) 2 (2g 12 - 49g n + 534g 10 - 3412g 9 + 14151g 8 - 39685g 7 

+75440g 6 - 92942g 5 + 62238g 4 - 455g 3 - 36019g 2 + 27008g - 6808) (8.7) 

A sqJ = 2q(q -2){q- 3)(q - l) 4 (g 9 - 19g 8 + 162g 7 - 818g 6 + 2706g 5 
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-6090g 4 + 9318g 3 - 9330g 2 + 5535g - 1477) 
A sq ,8 = -q(q - 4)(g - 2)\q - l) 5 (g 3 - 6q 2 + 13q - ll)x 

(g 4 - 9g 3 + 29q 2 - 40g + 22) 



(8.8) 



(8.9) 



B. Strip of the Kagome Lattice 

For our cyclic strip of the Kagome (kg) lattice we find 

A kgfi = q(q-l) 2 (q-2) 2 (8.10) 

A kgA = -q(q - l)(q - 2) 2 (2q 5 - 17g 4 + 59g 3 - 105g 2 + 104g - 59) (8.11) 

A kg ,2 = q(q - - 2)(2g 9 - 37g 8 + 311g 7 - 1547g 6 + 4994g 5 

-10845g 4 + 15932g 3 - 15484g 2 + 9250g - 2656) (8.12) 
A kg ,3 = q(q ~ 1)(<? - 2)(3q n - 69g 10 + 726g 9 - 4604g 8 + 19522g 7 - 58029g 6 

+ 123231g 5 - 186853g 4 + 198437g 3 - 141104g 2 + 60968g - 12292) (8.13) 
A kgA = -q(q - l) 2 (g - 2) 2 (2g n - 34g 10 + 251g 9 - 1020g 8 + 2290g 7 - 1840g 6 - 4218g 5 

+14163g 4 - 16324g 3 + 4750g 2 + 6204g - 4476) (8. 14) 
A kg , 5 = q(q - l) 2 (q - 2) 4 (7g 10 - 121g 9 + 953g 8 - 4521g 7 + 14396g 6 - 32498g 5 

+53581g 4 - 65213g 3 + 57416g 2 - 33432g + 9720) (8.15) 
A kgfi = -q(q - l) 3 (g - 2) 6 (9g 8 - 129g 7 + 810g 6 - 2927g 5 + 6708g 4 - 10133g 3 

+10142g 2 - 6446g + 2110) (8.16) 

A kgJ = q(q - l) 4 (g - 2) 8 (5g 6 - 54g 5 + 232g 4 - 511g 3 + 608g 2 - 380g + 118) (8.17) 

A kg , 8 = -q 2 (q-l) 6 ( q -2) u ( q -4) (8.18) 
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C. Asymmetric Homeomorphic Expansions of the Square Strip 

For the case k = 3 defined in the text, the coefficients entering in the numerator of the 



generating function are 

^=3,0 = q(q - 1)(9 - 2)(g 3 - 4g 2 + 7q - 5) (8.19) 

A k=3>1 = -q{q -l)(q- 2)(2g 4 - 12g 3 + 27q 2 - 29q + 13) (8.20) 

A k=3 ,2 = -q(q ~ lf{q ~ 2) V - q - 1) (8.21) 

^fc=3,3 = q(q ~ 2) 4 (g - l) 2 (g 2 - 2q + 2) (8.22) 
For the case k = 4 we find 

A k=4fi = q(q - l)(g 6 - 8g 5 + 28g 4 - 56g 3 + 70q 2 - 53q + 19) (8.23) 

A k=4il = q(q - l)(2q 7 - 21q 6 + 94g 5 - 236g 4 + 365g 3 - 356g 2 + 207g - 56) (8.24) 

A k=Ay2 = -q(q - l) 2 (2g 6 - I7q 5 + 59g 4 - 107g 3 + 107g 2 - 56q + 11) (8.25) 

A fc=4i3 = ~q(q -2)(q- l) 2 (q 3 - 5q 2 + 8q - 5)(g 4 - 5g 3 + 10g 2 - lOq + 5) (8.26) 
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